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ABSTRACT 


Non-thermal emission occurs in the cores of the 9.4 and 10.4Mm COg bands 
on Mars, and has been recently identified as a natural atmospheric laser. 

This paper presents observations of the total flux a*nd center-to-1 imb 
dependence of this emission for Mars and Venus. The emission is believed to 

be excited by absorption of solar flux in the near-IR COg bands, followed by 

collisional transfer to the 00°1 state of COg. A comparison is made between 
the observations and a detailed theoretical model based on this mechanism. It 
is found that the theoretical model successfully reproduces the observed 
center-to-1 imb dependence of this emission, to within the limits imposed by 
the spatial resolution of the observations. A comparison is also made between 
the observed fluxes and the predictions of the theoretical models. The 
observed flux from Mars agrees closely with the prediction of the model; the 

flux observed from Venus is 74% of the flux predicted by the model. 

We utilize this emission to obtain the kinetic temperatures of the 
Martian and Venusian mesospheres. For Mars near 70 km altitude, a rotational 
temperature analysis using five lines gives T » 135 ± 20K. The frequency 
width of the emission is also analyzed to derive a temperature of 126± 6K. In 
the case of the Venusian mesosphere near 109 km, the frequency width of the 
emission gives T * 204±10K. 
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I. INTRODUCTION 


Strong non- thermal emission in the 9.4 and 10.4Min bands, of COg was first 
reported for the atmospheres of Mars and Venus by Betz et al . (1976), Johnson 
et al . (1976) and Betz (1976) using infrared heterodyne spectroscopy. The 
emission occurs in the cores of the COp absorption lines and is only 

u g 

detectable at spectral resolving powers of >10 . The emission is believed to 
be formed by absorption of near- infrared solar flux, followed by colli si onal 
transfer of the absorbed quanta to the 00°1 level of COg* and radiative decay 
at lOum (Johnson et al . 1976; Betz 1976). The emission disappears in the 
absence of sunlight (Doming et al . 1982). Recently, Mumma et al . (1981) have 
demonstrated that in the case of Mars this emission arises from a population 
inversion, making it a natural laser. The emission intensity was reported by 
Betz (1976) and Betz et al . (1976) to be approximately proportional to tne 
incident solar flux, but little attempt was made to quantify this dependence. 
The understanding of this emission and its implications is still not complete* 
Until now, no complete comparison has been made between observations and the 
predictions of a suitably detailed theoretical model. 

This paper presents extensive observations of this emission from the 
mesospheres of Mars and Venus. We obtain the total flux emitted in the 9.4 
and 10.4}im bands and the center- to-limb dependence of the emergent intensity. 
These quantities are compared to the predictions of a detailed theoretical 
model (Deming and Mumma, 1983). We obtain the kinetic temperature in the 
Martian mesosphere near 70 km from the frequency width of this emission and 
from a rotational temperature analysis. 

II. OBSERVATIONS 

Observations of the emission in lines of the [(10°0)-(02°0)*J, I and II 
bands were made using the Goddard Space Flight Center infrared heterodyne 
spectrometer and the Kitt Peak National Observatory McMath main solar 
telescope during several observing runs from December 1979 through April 1980. 
The spectrometer is described in detail by Kostiuk et al . (1980) and Mumma et 
al . (1982). Observations were made of Mars near opposition and of Venus near 
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phase 0.5; the Venus observations (Fig. 1) are of higher quality than the Mars 
observations because the emission was much brighter and the planetary disk was 
much larger than the 1.7 arc-sec FWHM instrumental field of view, allowing 
better spatial resolution. Most of the observations were made in the 10.3337 
urn R(8) line, but other lines were also observed. The observations were made 
at two spectral resolutions: 25 MHZ (0.0008 cm"^) and 5 MHZ (0.000017 cm"^) , 
using two separate RF filter banks. For comparison, the Doppler half-width at 
half-maximum for COg lines at 150K is 19 MHZ. The 25 MHZ and 5 MHZ data were 
registered simultaneously; the 5 MHZ filter bank was tunable and was centered 
on the emission core. For all of these observations, terrestrial CO2 
absorption was present at the frequency of the COg laser local oscillator, and 
was defined and removed by observing the Moon, which is a thermal continuum 
source at this frequency. The observations were calibrated in terms of 
absolute flux, above the atmosphere, by adopting an appropriate temperature 
for the observed location on the lunar surface (Montgomery et al . 1966). For 
the planetary observations, the position of the beam was established by visual 
inspection of the image on a television monitor. Under conditions of good 
seeing we found that the beam center position can be set in this manner to a 
precision of «/«l arc-sec. This visually estimated position was later corrected 
for the (small) differential refraction between the visual and lOMm regions. 
This correction was made based on the indices of refraction for air at 
continuum frequencies. The aspect geometries of the planetary disks were used 
to compute the angular distance of each observed point from the sub-Earth and 
sub-solar points. 

III. ANALYSIS 

Each observation of the emission line yields three quantities: the 

frequency width of the line, tne line center frequency, and tne specific 

-2 -1 -1 

intensity of the emission core (ergs cm sec sr ). We obtain these 
quantities by fitting a gaussian to'the observed emission core, after 
stripping off the underlying absorption profile. This removal of the 
absorption profile was facilitated in the case of Mars by obtaining a model 
atmosphere fit to the absorption line. The absorption line profile was 
modelled as the sum of transmitted thermal radiation from the surface and the 
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self-emission of the overlying atmosphere* The surface temperature was 
obtained from the observed line profile by fitting to the flux In the far 
wings, where the atmosphere is optically thin. A temperature vs. altitude 
profile, appropriate to the local time and season, was adopted based on the 
discussion of Seiff (1978). The COg line parameters used in calculating the 
line profile are well known, and were taken from McClatchey et al . (1973) and 
Freed et al . (1980). The atmospheric surface pressure was treated as a free 
parameter and was varied in order to obtain an optimal fit to the observed 
line. Surface pressures derived in this way were consistent with the 
variations given in Viking ground-truth results (Hess et al. 1980). Examples 
of the modelling and removal of the underlying absorption profile are given on 
Fig. 2. Further discussion of the absorption line fitting process is not 
given here, because the results for the emission core are not sensitive to the 
details of this procedure. In the case of Venus the emission completely 
dominates the underlying absorption profile, which is so shallow that it can 
be ignored. 

a) Kinetic temperatures 


We derive kinetic temperatures for the emitting region by attributing 
the frequency width of the observed line to molecular thermal motion. In 
actual practice the derived kinetic temperatures are only upper limits, since 
other effects act to broaden the line. These effects include planetary 
rotation and atmospheric turbulence within the heterodyne field of view, 
possible local oscillator drift, and the finite resolution of the RF filters. 
During each observation there is also a change in the line of sight component 
of the earth's rotational velocity. The frequency widths derived for the 
emission core were corrected for this latter effect when it was appreciable. 

It is assumed in this procedure that gain-narrowing of the line profile due to 
stimulated emission is a negligible effect. The theoretical model of the 
emission (Doming and Numma 1983) supports this assumption. 

Kinetic temperatures were derived in this way only for lines in the 
10.4um band. For the better Mars data the error associated with a single such 
temperature determination is ±8K; the Venus data has smaller errors. However, 
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we find variations of »/»±30K from our Mars measurements, so we may be observing 
real temperature fluctuations (e.g. Zurek, 1976). Using this technique, we 
obtain average temperatures for the mesospheres of Mars and Venus of 126 ± 6K 
and 204 ± lOK at 70 and 109 km respectively. The individual measurements are 
given in Tables I and II. 

b) Rotational temperatures 

Since we expect the rotational levels of CO2 to be in thermal 
equilibrium at the pressures of relevance here, we can use observations of the 
emission from a range of rotational levels to determine the rotational 
temperature of the mesospheres. We have done this in the case of Mars by 
observing the emission from 6 lines with J-values ranging from 4 to 26. In 
the case of the R(26) line the emission core is very weak and the value we 
obtain for the integrated emission is sensitive to the manner in which the 
underlying absorption profile is modelled. This line was therefore omitted 
from the rotational temperature analysis, although if it were included it 
would- not -change the results significantly. ’ ’* ■ ~ . 

Under conditions of rotational thermal equilibrium the intensity of 
a rotational line is given (e.g. Herzberg 1950) as: 


Ij‘ ^ vSj' EXP (-B'J'(J'+1) hc/kT) 


( 1 ) 


where J‘ is the rotational quantum number for the upper state, v is the line 
frequency, and Sj^J' (R branch) or Sj=J'+l (P branch). We take B'*0.39 cm“^ 
(Freed et al . 1980). We normalize the observations to the intensity of the 
R(8) line having J* * 9. We use absolute Intensities for each line (discussed 
below) and we correct each intensity to its value at the sub-solar point using 
a theoretical model (Deming and Mumma 1983). The fluxes used in the analysis 
are given in Table I. The resulting rotational temperature for the Martian 
mesosphere is T^^^ » 135±20K, and the analysis is illustrated on Fig. 3. 
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c) Absolute flux measurements POOR QUAUTY 

It is important to derive an observed value for the absolute total 

flux emitted in the 9.4um and 10.4um bands. Theoretical modelling of the 

emission process can predict a value for thl$ flux, and comparison between 

observed and theoretical fluxes can shed light on the completeness of our 

understanding. In principle, the observed total flux can be obtained by 

observing the emission in each rotational line of each band and summing the 

observed values. This procedure is impractical because of the length of time 

required to observe a single rotational line with appropriate calibration. We 

therefore assume rotational thermal equilibrium at T»126K, as in (1), and we 

derive the total emitted flux by scaling that which is observed for a single 

line. We assume that the 9.4 and 10.4|im bands emit equally. This assumption 

is crudely consistent with our observations, and is prescribed more precisely 

by laboratory measurements of transition rates (Murray, Kruger and Mitchner 

1974). We note that the observations of each line give a specific intensity 

(ergs cm"'^ sec""^ sr”"^)', which we wish to integrate over solid angle to obtain 

-2 -1 

an emergent flux (ergs cm sec ). To accomplish this, we let * the arfgle 
between an observed ray and a line normal to the atmosphere. We let = the 
zenith angle of the sun at the observed position. We denote the emergent 
specific intensity as and we assume that (®£»®s) is independent 

of azimuthal angle An increment of solid angle, du, is given as du « 
sine^de^d* and so the emergent flux is 

^ ' J 

0 

We expect that I^( 9 £f 05 ) I^(O, 02 )/cos e^, 

hence we obtain 


»/2 

I^(0g,8s) sinogcos 0g d0g (2) 

0 



Vs) « (°’S) 


(3) 
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Observed values of given in Tables I and IL, Since we expect 

that the emission intensity will depend on solar zenith distance (e^), we use 
the predictions of the theoretical models to correct the observed valuas of 
l^(0g,es) to the sub-solar point {q^-O). We similarly correct these values to 
e£*0, and we give the resultant 1^(0, 0) values in Tables I and II. 

In the case of Mars we observed near opposition, and so ^ 0 at disk 

center. This means that the corrections to convert Iy( 0 £» 05 ) to 1^(0, 0) are 

small and the resultant 1^(0, 0) values are not significantly model dependent. 

For the 10.33um R(8) line we obtain 1^(0, 0) » 2.84 X 10"^ ± 0.34 X 10"^ ergs 

cm"*^ sec’^ sr“^. This gives f^(0) = 1.78 X 10’^ ergs cm"^ sec“^ for this 

line. We denote the total flux emergent in both bands as F (9^)' Adopting T « 

-2-1 vs 

126K, we calculate F^(0) * 13.7±1.6 ergs cm sec . Our observations, 
however, were made near aphelion (1.66 a.u.) and tha emission strength should 
be proportional to the incident solar flux. Correcting the emergent flux to 
the mean distance of 1.52 a^u. we obtain F^(0) - 16.2±i.9 ergs cm"" sec“*. 
Random errors in the observed emission intensities are the dominant source of 
uncertainty in F^(0). * 

In the case of Venus, the observing geometry creates more difficulty in 

correcting I^ (e^, o^) to I^ (0,0). In this case, it seems preferable to 

use only the 10.33um R(8) observations taken at 30° and 45° West. With this 

restriction we have (0,0) = 8.33 X 10 ergs cm sec sr . Adopting T = 

204K, from the kinetic width of the emission cores, we obtain F (0) = 56.2 
-2 -1 -1 v ' ' 

ergs cm sec sr . If, in contrast, we use all of the I (0,0) values for 

10.33nm R(8) we obtain I^ (0,0) * 6.93 ± 0.48 X 10“^ and F^(0) « 46.7 ± 3.2 

ergs cm”^ sec”^. The former value is preferable itince it will be less 

model -dependent and less prone to systematic error. 

d) Center-to-limb dependence of the emission 

One of the primary motivations for these observations was the hope that 
the center-to-limb dependence of the emission could be determined with 
sufficient accuracy to place meaningful constraints on theoretical models of 
the emission process. Our observations of the center-to-limb dependence are 


8 


given in Tables I and II as ratios of the emission intensity observed at each 
point relative to some normalization point • In the case of Mars this 
normalization point was taken at disk center. In the case of Venus the 
normalization point for the 10.33um observations was the 30° west observation. 
The 9.34um Venus observations were normalized to the average of the two west 
limb observations. For each observed point, Tables‘1 and II give values of cos 
Og and cos 0 g. We have also calculated theoretical values for the intensity 
at each point, relative to the normalization point. In calculating these 
theoretical ratios, Mars was represented by the 120K model. If we define I^ 
(0,0) » 1.0, inspection of the theoretical models showed that I^( 0 £i 05 ) « 

(0,Og) I^ ( 0 £»O). Values of I^ ( 0 , 02 ) and I^( 0 g,O) were tabulated in the 
theoretical models. Using this procedure we calculated I^Ceg,©^) for each 
point of the planetary disk which fell within the heterodyne field of view. 

The resulting theoretical intensities were convolved with the Airy pattern 
which represents the diffraction-limited beam pattern of the heterodyne 
spectroffleter. Modelled relative Intensities calculated in this manner are 
given in Tables I and il. 


The observed relative intensities are plotted versus the modeHed values 
on Fig. 4. The typical error bars which are given account for error in the 
observed relative intensities as well as error in the modelled values. 

Estimate of the latter error is based on a ± 1 arc-sec uncertainty in the 
position which was observed. Given the estimated random errors, the agreement 
between the observations and the model is good. However, systematic errors 
are more difficult to assess. Systematic errors could arise because we 
determined the position of the instrument beam by visually positioning the 
planetary image on a cross line reticle. Since subjective judgement is 
involved in this process, and since a variety of observers participated, 
systematic departure from our assumed ^trsitions is possible. 

Fig. 4 represents an improvement on the less quantitative comparisons 
given by Betz (1976) and Johnson et al. (1976). The principal value of Fig. 4 
is that it implies that the strength of the emission is indeed closely 
proportional to incident solar flux, i.e. that the theoretical values of I (0, 
0 ^) are correct. If this were not so, it would not have been possible to 
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obtain good agreement between the observed and theoretical Intensity ratios. 
With respect to the theoretical modelling of I^(0£,O) values, however. Fig. 4 
Is less useful. In particular, the observations do not have sufficient 
angular resolution to resolve the sharp peak in I (0p»O) which Is predicted to 
occur for cos e^^O.S. In this respect, a more significant test of the models 
could be obtained using data with much higher spatial resolution. 

IV. DISCUSSION 

Our measurements of the frequency width of the emission give a mean 
temperature of 126±6K for the Martian mesosphere near 70 km. This 
determination Is in agreement with our rotational temperature analysis, which 
yields 135±20K. These values are somewhat lower than the mean Viking probe 
results (Seiff and Kirk 1977), whlcn give ^139K, and the stellar occultation 
results of Elliot et al. (1977), which give 145±10K. In the case of Venus, 
our measurements of the thermal width of the emission give T » 204±10K near 
109 km, in agreement with the result from the Pioneer Venus Day Probe (Seiff 
' et al . 1980) at this altitude. " ..... 

-2 -1 

We derive a value of 16-2±1.9 ergs cm sec emitted at the sub-solar 

point in the 9.4 and 10.4iim bands from the Martian mesosphere. Our 

2 -1 

measurements for Venus give a total flux of 56±4 ergs cm sec emitted in 
both bands. These measurements are compared to the results of the theoretical 
models in Table III. The results of Johnson et al . (1976) for Mars are also 
included in this table. We note that Johnson et al . (1976) give a greater 
total flux than we obtain in this investigation, and they also obtained a 
higher kinetic temperature from the frequency width of the emission core. 
Johnson et al . (1976) derive T./*170K for the temperature near 75 km, a value 
which is significantly higher than other measurements have given. The 
theoretical models, however, show that their measured flux is quite consistent 
with their measured temperature. Unless this is a coincidence, it suggests 
that their determination of a 170K mesospheric temperature may represent a 
genuine mesospheric temperature fluctuation. 

The most significant conclusion which can be drawn from Table III is 
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that the emission observed from Mars is unexpectedly bright. The observed 
total flux is essentially equal, within the errors, to the flux predicted by 
the theoretical model. This is surprising, since Deming and Mumma (1983) note 
that the model makes assumptions which tend to overestimate the emitted flux. 
In the case of Venus, for example, the observed flux is only 74% of the 
predicted fVjx. Betz (1976) and Johnson et al . (1976) have suggested that 
near-IR absorption by water vapor, followed by resonant vibrational transfer 
to C02» can contribute to the pumping of this emission. This process is not 
included in the theoretical models, and may account for the differences shown 
by Table III. A final understanding of the 10^m CO2 emission is therefore 
tied tp a more comprehensive picture of chemical and radiative processes in 
the- mesospheres of Mars and Venus. 
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Figure Captions 


Fig. 1 

Fig. 2 

Fig. 3 
Fig. 4 


Heterodyne observations, at 5 MHz resolution, of emission in the 
10.33mm R(8) line of COg on Venus. The phase of the planet at the 
time of observation was 0.5 and the sub-solar point occurs at the 
(leftmost) limb. The intensity scale is normalized so that the lOum 
continuum (^235K) has a value of unity when observed at norm’.l 
incidence. The spatial resolution is indicated by the full width of 
the instrumental beam, to its half-power points (HPBW). 

Example shov/ing observations and modelling of the 10.33nm R(8) line 
of at the center of the Martian disk. The intensity scale 

is normalized so that unity .represents a brightness temperature of 
260K. The top portion shows 25 MHz data and modelled profiles; the 
bottom portion of the figure includes the 5 MHz observations of the 
emission core. 

Rotational temperature analysis of the laser emission from the 
Martian mesosphere. 

Comparison of theoretical and observational intensities for the 
laser emission. The data plotted is from Tables 1 and II. The 
abscissa gives the logarithm of the ratio of the intensity at an 
observed point to the intensity at a normalization point. The 
ordinate gives the theoretical value of this ratio calculated from 
the models. Typical error bars are shown; errors for individual 
points will vary depending on the details of the observing geometry. 
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TABLE I. Suiwary of Observations of the laser emission from Mars. Specific intensities are corrected for atmospheric absorption. The 

specific Intensity is a function of angular distance from the sub>earth point* and from the sub-solar point. e^. Nuobers in 
parentheses are powers of ten. 
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Table II« Suonary of Observations of the laser enission from Venus* Specific intensities are corrected for atmospheric absorption* The 
specific intensity is a function of angular distance from the sub-earth pointiQ c* >od from the sub-solar point. d ^* Numbers 
in parentheses are potters of ten. 
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Table 

III. Theoretical and observed 
9.4Mm and 10.4um bands. 

values of the total flux emergent in the 


Observed 

Theoretical 


ergs cm sec 

ergs cm**^ sec"^ 

Mars 

16±2 (X) 

15.1 (126K) (3) 


20 (2) 

20.6 (170K) 

Venus 

56±4 (1) 

75.3 


(1) This work 

(21 Johnson et al . (19761 

(3) Obtained by interpolating between the theoretical values at 120K and 
170K. 
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